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The analysis of intermittent data is improved. The standard method of 
recovering the history of a particle cascade is proved in general not to repro- 
duce the structure of the true cascade. The recovering corrections to the 
standard method are proposed and tested in the framework of multiplicative 
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cascading models. 



1 Introduction 



The first data on possible intermittent behaviour in multiparticle production [|TJ 
came from the analysis of the single event of high multiplicity recorded by the 
JACEE collaboration ||. Data from some accelerator experiments || confirmed 
afterwards that there are large dynamical fluctuations appearing in the high energy 
multiparticle final states which manifest a scaling behaviour. Many different models 
have been proposed since to explain the effect [|J. Some of them suggested that 
an underlying final state multiparticle cascade may be responsible for the scaling of 
multiparticle moments In this approach the intermittent data represent the last 
stage of the cascade, and the main problem lies in the extraction of the information 
on previous cascading stages which is in some way encoded in the last stage data. 
It should be stressed that the problem of recovering the history of the cascade may 
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not be solvable if adressed generally. However, the self-similar processes which are 
assumed to underly the final state structure obey some scaling law. This makes 
many features of local dynamics dissapear, and one may expect to extract from the 
last stage data at least a part of information on the real cascade parameters . 

The method of recovering the history of the cascade from the last stage data was 
proposed and applied originally to the JACEE event data. Since that time it became 
a standard tool of multiparticle data analysis @], especially in the event-by-event 
analysis M. 

In this paper we would like to improve the standard method of analizing the 
intermittent data, taking into account corrections due to the recovering the history 
of the particle cascade in the framework of multiplicative random cascading models 
|| [jj]. This problem has been already adressed and analized in part in Ref. ||. Our 
discussion will proceed as follows. In section 2 we characterize briefly the standard 
method of estimation of intermittency exponents, and introduce the definition of 
recovering corrections. In section 3 the definition of multiplicative models is pre- 
sented, and the special cases of multiplicative models: a—, p—, (p + a)— models 
are summarized. In section 4 the recursive equation for recovering corrections in 
a multiparticle model with possible neighbour-node memory is derived. Section 
5 is devoted to the implementation of recovering corrections into the analysis of 
data. The implementation algorithm is proposed and numerically tested. Finally, 
in section 6 we present our conclusions. 

2 Standard estimation of intermittency exponents 
and recovering corrections. 

Consider a sample of M bins describing an individual (intermittent) event. For 
simplicity assume that M = 2 n , where n is a natural number. We thus have 2 n 
numbers describing the content of each bin : 

xt\ « = 0,l,... l 2"-l (1) 

which represent e. g. the distribution of particle density into bins. One assumes that 
the bin ensemble has been generated in some cascading process, and the unnor- 
malized density moments for this process : 

manifest a scaling behaviour parametrized by intermittency exponents <p q : 

z (n) _ 2 n.<^ (3) 
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The standard method of estimation of intermittency exponents was introduced 
firstly for the analysis of JACEE event |2|. The method recovered the history of the 
cascade in the following manner: it established the value of density moments for each 
cascade step, and made the linear x 2 — fit to the points (k,logz^) (k = 1, 



,n, 



log x = log 2 x) : 

log + (4) 

In this way the eventual long-range correlations could be separated and would not 
contribute to the estimated slope <fi q . For the assumed bin-into-two-bins splitting 
scheme the true value of the x- n_fc ^ bin content was replaced by y| n : 

^ 2 k -l 

y% = ~^k~ x 2 k xi+j- ( 5 ) 

Z 3=0 

The intermittency exponents were extracted from the reconstructed moments 



(fc) 

*g; rec. 



2 k -l 

(k) 



Z q;rec. r>jfe 

Z 3=0 

assuming their power law behaviour : 

z {k) „ 2 k <. (7) 

q; rec. \ I 

So far (see e. g. |T^]) one has estimated the value of (normalized) <f> q , assuming 

Simply that (p q;n orm. = <P' q (<f>q;narm. ■= 4>q ~ 

There is, however, an open question if and how the cascade recovered from data 
refers to the true cascade which generated the data. For the purpose of estimating 
the intermittency exponents it is enough to ask about the relation between the true 
density moments z q and the reconstructed ones obtained from (P). It is obvious that 
formula (|5]) looses a piece of information on the primary cascade. We give a simple 
example to illustrate the problem. If we assume that the underlying cascading 
process preserves e. g. the total particle density Y^Jo~ x j = 1 (k = l,---,ri), then 
it follows from (|[) that the reconstructed cascade will not manifest this property: 
X^o 1 ^ 1- ft m eans that standard method does not recover the conservation law 
present in the true cascade. 

Moreover, it was found explicitly for the a— model H |7J (which does not preserve 
the total particle density) that there exists always discrepancy between the true value 
of 4> q and its estimation 0' (|7|), due to recovering technique (|5[). This problem was 
discussed in detail in Ref. f|TTf| . 

Assume we may express the discrepancy between the true and the reconstructed 
moments at the (n — k)th cascade step in a following way : 
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where the factor p g (k) denotes the corrections due to recovering procedure ([5]). We 
will call them recovering corrections. Corrections p q (k) contain information on the 
parameters of a specific process which generated the true cascade. It is obvious that 
they depend also on the cascade step. Substituting into (|) one arrives at the 
relation : 

log z%;*2 ~ Mp#)) — (n — k) ■ <J) q . corr . + b, (9) 

where the fitted slope (f) q]C orr. estimates the true intermittency exponent (fi q . 

In this paper we confine ourselves to recovering corrections considered for the 
class of multiplicative random cascading models. For the multiplicative cascade 
("multiplicative" means that at each cascade step the bin content is multiplied by 
a number to generate the bin content at the next cascade step) relation (||) holds 
explicitly, and recovering corrections take the form [[] l] (for proof see Appendix A) : 

p^) = (f^ 2 E4 fc) l ) (io) 

where the average (. . .) is taken over the (eventual) random choices while generating 
the cascade. The starting bin x^ is set equal 1. It is worth noticing that p q {k) may 
be also expressed in terms of the erraticity moments |§ : 

P q (k) = C 1 , q . (11) 



3 Multiplicative models 

As already mentioned, in our paper we restrict ourselves to the class of multiplicative 
random cascading processes with possible neighbour-node memory, generating the 
uniform distribution of particle density. The commonly used models of random 



cascading: a— model |12| and p— model |T3| belong to this class. For the purpose 



of testing our predictions for recovering corrections we introduce below a new many 
parameter (p + a)— model. 

In the multiplicative random cascading processes with possible neighbour-node 
memory assume for simplicity the root of a cascade to be equal to 1 : x^ = 1. One 
generates the next stages of the cascade recursively. The scheme is following. The 
two bins x^ 1 ^ and a^i+i are obtained from x\ k ^ by multiplication : 

:= W,.xf\ 

:= W 2 -xf\ (12) 
where W\ and W2 are random variables of the m model parameters aj, j = 1, . . . , m : 
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W\ = dj with probability p aj , 
W 2 = a,j with probability p aj , 

(13) 

with normalized probability weights p a . : 

m 

E^ = 1 - ( 14 ) 

The distribution of particle density will be uniform if the following condition is 
fulfilled : 

p(W 1 = ai, W 2 = a/) = p(W 1 = dj, W 2 = ai), (15) 

where p{Wi = a^W^ = %) denotes probability of choosing in flT2|) = a, and 
W 2 = cij = l,...,m). Then the density moments fulfill relation (|]), where 
intermittency exponents <p q are equal to : 

(f) q = log(a q lPai + ... + a q rn p am ). (16) 

The models : a—, p— and (p + a) may be derived from general multiplicative rule 
(|T2D. To obtain the a— model it is enough to assume random variables W±, W 2 to 
be independent : 

(W 1 W 2 ) = (W 1 )(W 2 ). (17) 

The a— model has no node memory therefore no conservation law can be imple- 
mented here. 

Relation (|TI|) reduces to the p— model after turning m = 2 and : 

a 2 = 1 — ai, 

P ai = Pa 2 = 0.5, 

p(W 2 = a 2 | Wi = ai) = 1, 

p (W 2 = a x | W t = a 2 ) = 1, (18) 

where p(W 2 = a, | W\ = aj) denotes conditional probability of W 2 = ai, if W\ = cij. 
The p— model is an example of a multiplicative model with the neighbour-node 
memory. Relation ( |18D implies that the sum x^ 1 ^ + ^2i+i = > an d the particle 
density in a node is preserved. 

Finally we introduce the many parameter (p + a)— model, using relation (|12j) 
combined with the a— and p— model restrictions : 

a 2i — 1 — a 2i-\i 
Van Pa,2i-\i 

(19) 
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where m is an even number (z = 1, . . . , ~), and : 

p{W 2 = a 2i | Wi = aat-i) = 1, 

p(W 2 = a 2i _x | W^i = a 2i ) = 1. (20) 

One may check that the particle distribution generated in the (p + a)— model is 
uniform. The (p + a)— model may involve any number of parameters Oj. Therefore 
it describes a more realistic case of cascading since for large m the distributions of 
particle density for W\ and W 2 fll3f ) may be approximated by a continuous distribu- 
tion function f{x) : W\ t2 = x with probability f(x) dx. The total particle density 
will be preserved for any m, according to (p0|). 



4 Recovering corrections in multiplicative models 

Now we calculate explicitly the correction p q (k) for any multiplicative model with 
possible neighbour-node memory. To do this we will split the bins (xj's) appearing 
in fllDD into a left half (i < 2 k - 1 ) and a right half (i > 2 k ~ 1 ) (T| : 



Using the fact that the left and right bins are independent, one arrives at the recur- 
rence equation: 

Pi( k ) = \ q E ( ) ) V 3 (k - l)p^(A; - l)(Wtwr j ) (22) 

which may be solved recursively together with the initial data : 

P„(0) = 1, 

Po(Ar) = 1. (23) 

A similar recurrence relation has been also obtained in ||. It should be stressed 
that coefficients (WlW^ 3 ) are the only parameters of the model needed to solve 
(|22|) recursively. It means that to calculate p q (k) for a given model we need only 
to know the coefficients (WlW^' 1 ) . In the next section we show how to apply this 
observation to the data analysis. 
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Finally we present recovering corrections calculated for the a— and p— models : 



P a q - m ° del (k) = ^ £ (J J Pi(k - l)p q -j(k ~ 1)2*-*^ , (24) 



'q 

where <pj denote intermittency exponent ( fl6| ) 



5 Implementation of the recovering corrections 

The idea how to implement the corrections p q (k) into the analysis of the a— model 
data was sketched briefly in |Tl |. Here we extend the primary scheme, and apply 
it to the multiplicative model data. As it was mentioned in the previous section, 
coefficients (WlW^ 3 ) are the only parameters of the model needed to calculate 
recursively the corrections p q {k). Let us introduce a new notation for (WlW^) '■ 

{W{W l 2 ) = k hl . (26) 

We ask now how to derive kj q _j from the model. One may notice that for either 
j = or I = coefficients kjj equal : 

/.-,o k (lj ■>'■■. (27) 

where 0j's are ordinary intermittency exponents (§) which may be determined from 
relations (|J), ([|). To find the value of kjj ^ 0) we use the unnormalized 
density correlators off [[l], [|, |IIJ : 



c?J = -=- v UVVUV.Y 

3, 2 k . ( . 



In multiplicative models the correlators and the density moments fulfill the relation : 

of} = 47 1} ■ hi (29) 

which can be also rewritten as : 

hgcfj = (k-l)(f> j+ i + logk^. (30) 

Relation (|29"D may be easily derived from (|l|) : since each term in sum (]28|) originates 
from one node xf~ l \ it can be rewritten as : (x^Y (4*+i) = (^'^Y W(wl 
Relations (£2^), ( POj ) imply that we may derive k^\ in a straigthforward way by 
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calculating correlators and density moments from data, and applying to them the 
standard x 2— fit- 
Applying the standard method to the correlators at the previous cascade stages, 
we expect to find the discrepancy (due to the recovering procedure) between re- 
constructed correlators and the true ones, similarly as for the density moments. It 
can be proved (see Appendix B) that the discrepancy may be expressed in terms of 
recovering correction p q (k) fllOD : 



(n—k) (n—k) 



rec. ^j,l 



p j+l (k). (31) 



Now we have all tools needed for implementation of recovering corrections into the 
multiplicative data analysis. Below we propose an implementation algorithm which 
recursively adjusts the primary parameters <fi q , k^i (j + 1 — q, jl > 0) obtained after 
applying the standard method to the data : 

(INPUT) parameters cf)^ . . . , tf> q _i, kjj (j + I = 1, . . . , q — 1) obtained after 
applying the implementation algorithm for q = 1, 2, . . . , q — 1 step-by-step 
(for determination of 0i see Appendix C.3), 

(1) derive (f)' q , k'- q _ - {j — 1, . . . , q — 1) from data, using the standard method 
i. e. reconstruct the cascade using @ and derive the parameters from relations : 

log^fL =*•< + &. (32) 

C j,l;rec. ~ z j+l;rec. \°°J 

where k = l,...,n (cf. (|), ©), 

(2) derive (t> q - corr , , fcj, g -j ;C orr. (j = 1, • • • , q — 1) in the following substeps : 
(2.0) calculate p q (k) from relation (cf. (p2|)) : 

P«( k ) = J q E ( ) ) Pj(k ~ l)p 9 -j(k - l)k hq . v (34) 

using (f)' q , k'j^j derived in step (1), and estimate 4> q . corr , from (cf. @) : 

log 4-re k J. ~ MP#)) — (n — k) ■ <P q , corr . + b, (35) 
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(2.1 ) calculate p q (k) from using (p q - CO rr. (other parameters as after step (1)), 
and estimate corn from relation (cf. (|30|) , fl31~D , see also Appendix C.4) : 



.(»-*) 

rec. 



log(pj+j(fc)) = {n-k- l)4> j+ i + log fcj, j 



(36) 



(2.q-l) calculate from (|34]), using all previously derived parameters cf) q . corr , 

kj,q-j-,corr. > an d estimate feq-1,1; corr. from (|36"D, 

(3) compare the values of <f>' q , k' jq _ 5 and <f> q . corr ., k jiq _j. corr , (j = 1, . . . , q - 1). If 
the relative difference is large, assume : 



and repeat steps (2), (3) recursively until the relative difference between parameters 



(OUTPUT) parameters </> 1; . . . , g ; fe^j (j + 1 = 1, . . . , q). 

Other techniqual details and problems which may appear when applying the al- 
gorithm to data are listed in Appendix C. 

We have performed numerical simulations of the a—, p— and (p + a)— models 
in order to test how the implementation algorithm works in practice. We generated 
10000 cascades of the 10 step length for the a— and (p+a)— models, and one cascade 
of the 10 step length for the p— model [] for two different parameter sets separately. 

Implementation algorithm analized the data of the last cascade step. For each 
event it estimated the value of normalized intermittency exponents 4>2;norm., (p3;norm. 
((pi; norm. '■= 0i — using the standard method (step 1) with recovering corrections 
included (steps 2,3). The results are presented in Figs. 1, 2, 3, 4 (for the a— and 
(p + a)— models) and in Tabs. 1, 2. 

For the a— model the histograms of 2; norm ., <p3- n orm. obtained in the standard 
method and the histograms with recovering corrections included are almost identical. 
In this case the recovering corrections can be implemented better when one applies 
directly dedicated a— model recovering correction p^- model (k) (§D (see Figs. 1, 2 
and Tab. 1). The different accuracy of these both approaches is due to the fact that 

it can be proved that for a given parameter set the model generates always the same values 
of the correlators and density moments 
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the random variables W±, W% are independent in the a— model, and the coefficients 
kj >q -j ( p6[) are approximated better by the product 2^> • 2 < ^ _J ' than from correlators 




On the contrary, the implementation algorithm works well for the (p + a)— model 
(see Figs. 3, 4 and Tab. 2). For the (p + a) — model the histogram with the recovering 
corrections included approximates well the theoretical value of normalized intermit- 
tency exponent. The histogram obtained by using the standard method is moved 
slightly to the left in comparison to the histogram with recovering corrections in- 
cluded. 

We have checked that for the p— model the theoretical values of normalized in- 
termittency exponents are estimated perfectly by both standard method and imple- 
mentation algorithm, as we have expected 0. 

It should be also mentioned that the histograms generated by both implemen- 
tation algorithm and dedicated recovering corrections are symmetric, in contrast to 
the standard ones, and their dispersions are relatively small (see Tabs. 1,2). Finally, 
we notice that the accuracy of the estimation of intermittency exponents, while ap- 
plying the standard method and the improved one to a given model, depends on the 
parameters of this model (cf. Figs. 1, 2 and Figs. 3, 4). In Appendix D we present 
a qualitative analysis of the effect for intermittency exponents of the second rank. 
A similar analysis has also been done in ||. 

6 Conclusions 

To sum up we analized the estimation of intermittency exponents from the data 
which were generated by a multiplicative random cascading process. The following 
methods were applied: the standard method of cascade recovering ([|) and the im- 
proved method which included recursively the recovering corrections. The improved 
method was applied in the form of the implementation algorithm. Numerical sim- 
ulations have been performed to check how both methods work in practice. The 
conclusions may be summarized as follows : 

(a) standard method of estimation of intermittency exponents does not apply 
for the whole class of multiplicative models : its accuracy depends on the specific 
properties of the model and its parameters. The method does not detect a conser- 
vation law if present in the model; 

(b) improved method of estimation of intermittency exponents applied in the 
form of recursive implementation algorithm either corrects the standard method 

2 it follows from relations (p|),(p5[) that in the p— model 4>i ]T ec. ,norm. — 4>i;norm. if i > 1 
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estimation or does not change the standard method result. In the latter case the 
estimation may be corrected by applying the dedicated recovering corrections. In 
any case the improved method tests the applicability of the standard method, and 
allows one to estimate the accuracy of the intermittency exponents estimation. 

While applying the improved method, the parameters of the model are adjusted 
recursively from the primary (standard method) parameters. The histograms gen- 
erated by the improved method are symmetric, and their dispersions are of the 
same order as those determined for the standard method. The improved method 
takes into account the neighbour-node memory (a conservation law) if present in 
the model, by including the density correlators into the estimation scheme. 
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Appendix A 

We prove relation (fL0|). The density moment z^ n ~ e ^) may be rewritten as : 



(n-k)_J_ v ( y {n - k A q = — V iVx (n) 1 (37) 

Z i=0 Z i=0 \ Z j=0 



Notice that : 

(n) _ (n-k) (k) 

X 2 k xi+j ~ X i ' 00 j ■ 

Substituting fl3"E| ) into (0), one arrives at : 

1 ^ / (n-kWI / 1 



>-*) 

"q;rec. nn-k ^ J \ Ok 

z i=o V z j=o 



Appendix B 

We prove relation (|31~D. The correlator c^"^ may be rewritten as 

^ IT {&-»)' '(&")' 



2" 

(n-k) _ _ 
rec. i^r, 

i=0 



(3f 



~ k e (*t k) y iE-f =*r°p,(*)- (39) 



ii 



2-1 \ / , 2-1 



-fe-1 ( 9fc X 2*x2i+m ) ( 9fe X 2 fc x (2i+l)+r J - (^) 

Relation (RSI) implies : 



2n-fe— 1 / i I 2 fc / < 2 K x2i+m I I 2« 

i—O V m— / \ r— 



cfci' = ^rr r E" (4r*>y (*£?>)' (? = <&-*W*>- (4i) 

z t=o V z i=o / 



Appendix C 

We list some techniqual details which can be useful when applying the implementa- 
tion algorithm to data, and discuss possible problems. 

(1) The recovering corrections applied for calculating 2 , 03 read : 

Pl (k) = 2 fc01 (42) 
P2 (k) = -{p 2 {k-l)k 2 , +p 2 {k-l)k 0t2 + 2p 2 1 {k-l)k 1A ) (43) 

Pa(k) = ~{p3{k - 1)^3,0 + P3(k - l)feo,3 

+ 3p 1 (k-l)p 2 (k-l)h t2 + 3p 2 (k-l)p 1 (k-l)k 2<1 ) (44) 



2) Since random variables Wi, W% generate the uniform distribution (cf. dl5|)) 

(45) 



relations (^,^,0) may be simplified by substituting : 

kj,i = kl,j- 



The (experimental) estimation of kj t i will be better if we determine c^f as 



(n) 



i=0 



(46) 



(3) Determination of It follows from relations (p|),(|42j) that the unnormalized 
reconstructed density moment Zi- rec . for any multiplicative model takes the form : 

z[% c . =2^=z{ n \ (47) 

which does not depend on the cascade stage k. Therefore the proposed estimation 
of 0i, implied by fl47D reads : 

= logjgg, 
n 
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(4) Estimation of coefficients kjj in the implementation algorithm. To estimate 
coefficients kjj in the step (1) of the implementation algorithm we calculate the 
reconstructed correlators and density moments, and apply relation (|33|) whereas in 
the step (2) we use relation (|36| ) to do the same. 

The brief explanation of the method is following. In formula (|33"D kjj appears as 
a slope, and an ordinary linear x 2_ fit may estimate it with a good accuracy. This 
approach works well for the reconstructed correlators and moments. 

On the contrary, the long-range correlations: log b present in relation (^) add to 
the value of log k^i, and generate a large error while estimating kjj := exp(log kjj + 
log b) from (|36|) . We could try to estimate b, assuming that : 

z (n) = 2 n-<f> qb ( 49 ) 

Then it follows from (|36| ) that : 

log cj? = (n - l)(f>j+i + log k jh (50) 

and log kj t i equals : 

log k jt i = log kj t i - log b. (51) 

The latter approach does not work for the reconstructed moments, where relations 
(^5|,[5H|) apply only approximately. However, it applies quite well for the moments 



with recovering corrections included because for this case relation (|3^) would require 
including recovering corrections to both correlators and density moments which in 
turn would generate larger error in estimating kjj. 

It was checked that the above method works for the multiplicative model data. 
However, the problem of determination of coefficients kjj and, in particular, the 
problem of the determination of correlators c^" from the real data is much more 
complicated (see e. g. p, pp. In this case the method needs some improvement 
which we will not discuss here. 



(5) The recursive implementation algorithm is not always convergent. Since 
estimation at the qth step depends upon parameters which were adjusted in the 
previous steps l,...,q — 1, the estimation errors propagate and get larger with 
growing q. Then it happens sometimes (not very often), that recursive adjusting 
ends with the repeating a sequence of different values of parameters, or parameters 
become not definite. In such a case we stop the algorithm, assuming for the values 
of intermittency parameters those derived in step (1). 



Appendix D 

The accuracy of the estimation of intermittency exponents, while applying the stan- 
dard method and the improved one to a given model, depends on the parameters of 
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this model (cf. Figs. 1, 2 and Figs. 3, 4). Below we present a qualitative analysis of 
the effect for the intermittency exponents of the second rank. 

Equation fl43|) may be solved analytically, and the solution (valid for any multi- 
plicative model) takes the closed form : 

p 2 (k) = (1 - A) ■ 2 { ^- 1)k + A ■ 2 2fc ^ , (52) 

where : 

A = 0A jfe - . (53) 
The reconstructed moments then read : 



z. 



(k) 

2; rec. 



m _ 2^ 2 ~ 1 ^ n . 2 k + A2^ in ■ 2^ 2 ' ,7lorm - k (54) 



There are two power law terms : 2 k and 2^ 2 ' norm - k in z^\ ec , ■ In order to establish 
how they influence the determination of <p' 2 (§,0) we performed the following check. 
For a given multiplicative model with fixed parameters (e. g. a— model with param- 
eters as in Figs. 1, 2) we established the values of z$ rec , from (H), and made the 

linear \ 2 ~ fit to the points (k, log z^l ec , )■ The slope obtained from the fit estimated 
the value of 4>' 2 . For the a— model we have obtained : for case a) <p' 2 = 0.0250, and 
for case b) <p 2 = 0.284. Both results agree with the 0' 2 's obtained from the model 
simulation (cf. Figs. 1,2 and Tab. 1). Similar analysis can be also peformed for the 
(p + a)— model. 

The above results confirm our observation that the accuracy of the estimation of 
intermittency exponents from the standard method depends on the parameters of 
the model. And, if the estimation of primary parameters (p' q , kj t is more accurate, 
their recursive adjusting performed by the implementation algorithm will be faster 
and more accurate as well. It means that in this case the improved analysis works 
better as well. 
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Figure captions 

Figs. 1, 2 Estimation of normalized intermittency exponents <fi 2 - n0 rm. and 03 ;norm . 
for a— model, using the standard method (dotted line), the improved method with 
the implementation algorithm (thin solid line), and dedicated a- corrections 
(dashed line) compared with the theoretical values (solid line), performed for two 
sets of a— model parameters : 

a) a± = 0.8, a 2 = 1.1, p\ = 1/3 

b) ai = 0.5, a 2 = 1.5, pi = 1/2. 

Figs. 3, 4 Estimation of normalized intermittency exponents <fi 2 - n0 rm. and (ps-norm. 
for (p + a)— model, using the standard method (dotted line), the improved method 
with the implementation algorithm (thin solid line), compared with the theoretical 
values (solid line), performed for two sets of (p + a)— model parameters : 

a) a 2i = 1 - a 2i -i, p 2i = P2*-i = 0.05 for % = 1, . . . , 10, 
a\ = 0.2, ci3 = 0.5, a 5 = 0.6, a 7 = 0.3, a 9 = 0.45, 

an = 0.25, ai3 = 0.1, ai5 = 0.15, an = 0.87, dig = 0.66, 

b) a 2i = l- a 2 i-i, £>2i = P2i-i for i = 1, . . . , 10, 
a,\ = 0.2, a 3 = 0.5, a 5 = 0.6, a 7 = 0.3, a 9 = 0.45, 

a u = 0.25, 013 = 0.1, ai 5 = 0.15, 017 = 0.87, a 19 = 0.66, 

2 Pl = 0.05, 2j9 3 = 0.15, 2p 5 = 0.25, 2p 7 = 0.40, 2p 9 = 0.05, 

2pn = 0.05, 2pi 3 = 0.02, 2p 15 = 0.02, 2p 17 = 0.005, 2p 19 = 0.005, 
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Tab. 1 Estimation of normalized intermittency exponents 4>2-,norm. and 03 ;norm . and 
their dispersions for the a— model, using the standard method (second column), the 
improved method with the implementation algorithm (third column), and dedicated 
a— corrections (p4|) (fourth column), compared with the theoretical values (first 
column), performed for two sets of a— model parameters (see Figs. 1, 2) : 





theor. 


standard 


algorithm 


a-corr. 


a) 02; norm. 


0.0285 


0.0251 ±0.004 


0.0246 ± 0.0033 


0.0288 ± 0.004 


03; norm. 


0.0813 


0.0757 ±0.010 


0.0727 ±0.009 


0.0798 ±0.0111 


b) (f>2;norm. 


0.322 


0.264 ± 0.044 


0.253 ±0.050 


0.276 ±0.051 


03; norm. 


0.807 


0.653 ±0.105 


0.666 ±0.125 


0.750 ±0.131 



Tab. 2 Estimation of normalized intermittency exponents 4>2-norm. and 4>^ n orm. and 
their dispersions for the (p+a) — model, using the standard method (second column), 
the improved method with the implementation algorithm (third column), compared 
with the theoretical values (first column), performed for two sets of (p ± a)— model 
parameters (see Figs. 3, 4) : 





theor. 


standard 


algorithm 


a) 02; norm. 


0.333 


0.322 ± 0.044 


0.305 ±0.066 


03; norm. 


0.832 


0.736 ±0.118 


0.720 ±0.174 


h) 02; norm. 


0.177 


0.170 ±0.023 


0.173 ±0.029 


03; norm. 


0.478 


0.438 ± 0.069 


0.470 ±0.092 
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Figure 1: 
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Figure 2: 
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Figure 3: 
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Figure 4: 



